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Abstract 

PpcA is the most abundant member of a family of five triheme cytochromes Cy In the bacterium Geobacter sulfurreducens 
(Gs) and Is the most likely carrier of electrons destined for outer surface during respiration on solid metal oxides, a process 
that requires extracellular electron transfer. This cytochrome has the highest content of lysine residues (24%) among the 
family, and It was suggested to be involved in e~/H^ energy transduction processes. In the present work, we investigated 
the functional role of lysine residues strategically located in the vicinity of each heme group. Each lysine was replaced by 
glutamine or glutamic acid to evaluate the effects of a neutral or negatively charged residue In each position. The results 
showed that replacing Lys^ (located near heme IV), Lys^^ (near heme I) or Lys^^ (between hemes I and III) has essentially no 
effect on the redox properties of the heme groups and are probably Involved In redox partner recognition. On the other 
hand, Lys'*^ (near heme IV), Lys^^ (between hemes III and IV) and Lys*° (near heme III) are crucial in the regulation of the 
functional mechanism of PpcA, namely in the selection of microstates that allow the protein to establish preferential e~/H^ 
transfer pathways. The results showed that the preferred e~/H^ transfer pathways are only established when heme III is the 
last heme to oxidize, a feature reinforced by a higher difference between its reduction potential and that of its predecessor 
In the order of oxidation. We also showed that K43 and K52 mutants keep the mechanistic features of PpcA by establishing 
preferential e~/H^ transfer pathways at lower reduction potential values than the wild-type protein, a property that can 
enable rational design of Gs strains with optimized extracellular electron transfer capabilities. 
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Introduction 

Geohacler species are of special interest because of their 
capability of transferring electrons to extracellular substrates 
during anaerobic respiration on insoluble terminal electron 
acceptors. This feature is currently being explored in developing 
Geobacter-h?Lsed biotechnological applications, particularly in the 
fields of bioremediation and bioenergy, such as removal of toxic or 
radioactive metals and electricity production by microbial fuel cells 
[1-3]. The detailed comparative analysis of the genomes of six 
Geobacter species (G. bemidjiensis, G. strain FRC-32, G. lovleyi, G. 
metallireducens, G. sulfurreducens and G. uraniireducens) showed 
that, except for G. lovleyi, the most remarkable feature is the 
presence of a large number and diversity of c-type cytochromes 



[4-7]. However, most cytochromes are poorly conserved among 
the Geobacter species [4] . One noticeable exception to this general 
feature is the family of triheme periplasmic cytochromes Cy present 
in all species of Geobacter. The triheme cytochromes from G. 
sulfurreducens (fis) is the only family studied in detail to date by 
genetic, functional and structural methods [8-18]. This family 
consists of five cytochromes, designated PpcA-E with approxi- 
mately 10 kDa molecular weight [14]. Each protein contains three 
low-spin heme groups with bis-histidinyl axial coordination and a 
conserved overall fold [11,18]. The functional electron transfer 
mechanisms of Gs triheme cytochromes have been characterized 
in detail, except for PpcC due to the presence of multiple 
conformers during the redox cycle of this protein [19]. The results 
obtained suggested that PpcA and PpcD can couple e /H"*" 
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transfer, a property that might contribute to the gx-neration of a 
proton electrochemical gradient across the periplasmic membrane 
[17,20]. Interestingly, PpcA and PpcD are the only members of 
the family whose deletion has significantly altered Gs phenotype 
with respect to the reduction of iron oxides [15,16]. 

A salient feature of Gs triheme cytochromes is their high lysine 
content as shown in Figure S 1 . This feature is expected to have a 
role in the modulation of the functional properties of these 
cytochromes, as well as in the recognition of redox partners. 
Amongst this family, PpcA is the most abundant cytochrome in 
the Gs periplasm and has the highest content of lysine (24%). In 
addition to its crystal structure [14], the solution structure of PpcA 
was determined [1 1] and can be used for the rational design of 
PpcA mutants. In the present work, we aimed to probe the 
functional role of selected lysine residues in the vicinity of each 
heme group in PpcA. The spatial localization of the selected lysine 
residues in PpcA structure is depicted in Figure 1: lysine 18 (near 
heme I), lysine 22 (between hemes I and III), lysine 60 (near heme 
III) and lysines 9, 43 and 52 (near heme IV). Each lysine was 
substituted by glutamine and glutamic acid to evaluate the effects 
of replacing the positively charged amino acid at each position by 
a neutral or negatively charged residue. For each of the twelve 
mutants, the redox properties were first screened at pH 6 and 
pH 8. For those mutants showing clear differences in redox 
behavior compared to the wild-type cytochrome, the detailed 
thermodynamic properties and the functional electron transfer 
mechanisms were determined using visible and NMR spectrosco- 
py. The comparison of the redox properties of the mutants with 
those of the wild-type cytochrome allowed us to elucidate the role 
of the selected lysine residues in the functional mechanism of 
PpcA. 

Materials and Methods 

Site directed mutagenesis 

For mutagenesis, QuikChange Site-Directed Mutagenesis Kit 
(Stratagene) was used in accordance with the manufacturer's 
instructions. Oligonucleotides were synthesized by MWG Biotech 
(High Point, NC). PpcA expression vector pCK32 [21] was used as 

a template. The presence of desired mutations was confirmed by 
DNA sequencing performed by MWG Biotech. 

Bacterial growth and purification of PpcA mutants 

Escherichia coli strain BL21(DE3) containing the plasmid 
pEC86 [22] was transformed with the expression vector contain- 
ing each of the PpcA mutant sequences (K9Q/E, K18Q/E, 
K22Q/E, K43Q/E, K52Q/E and K60Q/E). Transformed E. 
coli cells were grown in 2xYT medium containing 34 |Xg/mL 
chloramphenicol and 100 |ig/mL ampiciUin. Protein expression 
was induced by adding isopropyl-fi-D-thiogalactoside (IPTG) to a 
final concentration of 10 |xM. PpcA mutants were purified by 
cation exchange and gel filtration, as described for the wild-t\pe 
cytochrome [14]. The purity of the proteins was evaluate by SDS- 
PAGE (15%), stained with Coomassie blue. To confirm the correct 
heme incorporation, electrospray ionization mass spectrometry 
was performed at the HHMI Biopolymer Laboratory and W.M. 
Keck Foundation Biotechnology Resource Laboratory at Yale 
University. 

NMR studies 

Preparation of NMR samples. NMR samples and exper- 
imental conditions matched those used in the characterization of 
the wild-type protein in the reduced and intermediate state of 
oxidation [17] and are summarized below. PpcA lysine mutant 



sampk-s of about 140 nM for fuUy reduced studies and 70 |a,M for 
redox titrations were prepared from lyophilized protein in 80 mM 
phosphate buffer with NaCl (250 mM final ionic strength) in D2O. 
Full reduction of the samples was achieved by adding gaseous 
hydrogen in the presence of catalytic amounts of the enzyme Fe- 
hydrogcnasc-, from Desulfovihrio vulgaris (Hildenborough). Par- 
tially oxidized samples were obtained by first removing the 
hydrogen from the reduced sample with argon followed by the 
addition of controlled amounts of air into the NMR tube. The pH 
of the samples at the reduced and intermediate stages of oxidation 
was adjusted inside an anaerobic glove chamber with argon 
circulation to avoid sample reoxidation. 

NMR experiments and assignment. All the NMR exper- 
iments were performed on a Bruker Avance 600 spectrometer 
equipped with a triple-resonance cryoprobe (TCI) at 288K. As in 
the case of the wild-type cytochrome [17], NMR spectra were 
obtained before and after the lyophUization of each mutant to 
confirm that the protein integrity was not affected. NMR 
experiments for the fuUy reduced samples were acquired in the 
same experimental conditions to those of wild-type protein 
[10,23]. The following experiments were acquired: 2D-'H 
TOCSY (45 ms) and 2D-'H-N0ESY (100 ms). To monitor the 
stepwise oxidation of the individual hemes in PpcA lysine mutants, 
2D-'H-exchange spectroscopy (EXSY) spectra in the pH range 6— 
9 were acquired in partisdly oxidized samples, as previously 
described for the wild-type protein [17]. Spectra were processed 
using TOPSPIN (BrukerBiospin, Karlsruhe, Germany). The 
program Sparky - NMR Assignment and Integration Software 
[24] was used for inspection and for assignment of the signals in 
the NMR spectra. The assignment of specific heme signals to the 
corresponding hemes in the reduced state of PpcA lysine mutants 
was performed using the methodology previously used for the 
wild-type protein [10,17]. 

Redox titrations followed by visible spectroscopy 

Due to the negative reduction potentials of the redox centres, 
PpcA lysine mutant redox titrations followed by visible spectros- 
copy, were carried out inside an anaerobic chamber (MBraun) as 
described previously [10]. Samples with 18 (xM protein concen- 
tration were prepared at pH 7 and 8 in NaCl/phosphate buffer 
solution with a final ionic strength of 250 mM. Redox titrations 
were carried out at 288 K. The reduced fraction of the proteins 
was determined by integrating the area of the a-peak (551 nm) 
above the line connecting the flanking isosbestic points (543 and 
559 nm) to subtract the optical contribution of the redox 
mediators, as described for the wild-type protein [10]. To check 
for hysteresis, each redox titration was performed in both oxidative 
and reductive directions. The experiments were performed at least 
two times, and the reduction potentials (relative to standard 
hydrogen electrode, SHE) were found to be reproducible within ± 
5 mV. 

Thermodynamic model 

In the case of a triheme cytochrome, three consecutive 
reversible steps of one-electron transfer convert the fuUy reduced 
state (stage 0, Sq) to the fully oxidized state (stage 3, S3), and 
therefore four different redox stag(;s can be defined, each group 
consisting of microstates with the same number of oxidized hemes 
(see Fig. S2), as previously described for the wild- type protein [17]. 
The general theoretical framework that allows the detailed study of 
the redox centre thermodynamic properties in multiheme proteins 
[25] was used to determine the thermodynamic parameters of the 
redox centres for the wild-type protein [17] and is summarised 
below for clarity. In summary, taking as reference the fully 
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Figure 1. Spatial location of lysine residues mutated in this work depicted in PpcA solution structure (PDB code, 2LDO [11]). Tine 
PpcA polypeptide chain (gray) is shown as Q ribbon and heme groups (red). The side chain of Lys', Lys'^ Lys'^^ Lys''^ Lys^^ and Lys^° (green) are 
represented as stick drawings. The hemes are numbered I, III and IV, a designation that derives from the superimposition of the hemes in 
cytochromes c? with those of the structurally homologous tetraheme cytochromes C3. 
doi:10.1371/journal.pone.0105566.g001 



reduced and protonated protein, the energy of the microstates can 
be described in the full range of pH and solution potential as sums 
of 10 parameters: the three energies of oxidation of the hemes 
(reduction potentials), the pK^ of the redox-Bohr centre, plus six 
two-centre interactions energies (three heme-heme and three 
redox-Bohr interactions). To obtain such information it is 
necessary to monitor the stepwise oxidation of each heme using 
2D- H EXSY NMR experiments, which allows discriminating the 
individual heme signals in different oxidation stages [26,27]. The 
information obtamed from the NMR data allows the determma- 



tion of the thermodynamic parameters relative to reference state 
(for a review see [25]). 

In the characterization of the wild-type protein the chemical 
shifts obtained for the heme methyls 12'CH3', l^CU-/^, 
12'CH3'^ (labelled according to the lUPAC nomenclature - see 
Figure S3) in each oxidation stage in the pH range 6-9 were fitted 
simultaneously with visible redox titrations data obtained at pH 7 
and 8. 

In order to obtain thermodynamic parameters of the PpcA 
lysine mutants the same set of heme methyls were used and visible 
redox titrations were performed in identical experimental condi- 
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tions. The experimental uncertainty of the NMR data was 
evaluated from the line width of each NMR signal at half height; 
the visible data points were given an uncertainty of 3% of the total 
optical signal. 

Results and Discussion 

Impact of the mutations on the cytochromes global fold 
and heme core 

The yields of the twelve PpcA mutants were similar to that 
obtained for the wild-type protein (approximately 3 mg per liter of 
culture). The mutants showed identical UV-visible and NMR 
spectroscopic features both in the reduced and oxidized states 
compared to the wild-type protein (data not shown). 

Due to the small size of triheme cytochrome PpcA (71 residues), 
the core of the protein formed by the three heme groups is 
essentially covered by the polypeptide chain in an approximate 
ratio of 24 residues per heme group. Moreover, due to the strong 
ring-current efferts produced by the heme groups, the proton 
chemical shifts of the heme substituents are determined by the 
relative orientation of the hemes and neighbouring residues. 
Therefore, heme proton chemical shifts provide an excellent probe 
to detect major changes in the overall fold of the polypeptide chain 
and to evaluate the impact of the substitution of each lysine by 
glutamine or glutamic acid. The impact of the mutations on the 
heme core architecture was probed by 2D- H NMR. The heme 
proton resonances were assigned using the same methodology 
described for PpcA [10] and their chemical shifts were compared 
to those of wild-type cytochrome (see Table SI and Fig. S4). The 
rmsd values calculated from these chemical shifts are low, though, 
as expected slightly higher for the heme signals that are closer to 
the mutated residue. The small values of rmsd and, therefore, the 
good correlation obtained for the chemical shifts of heme protons 
in the native and mutated cytochromes, indicate that the heme 
core of the proteins were unaffected by the mutation. In addition, 
the same set of NOE connectivities between the heme groups were 
observed for native and mutant proteins, further confirming that 
the heme core arrangement is similar in all the proteins. 

In silico modeling of mutations 

Side chains of lysine residues that are fully exposed to solvent 
are usually not well-defined and can have many different 
conformations resulting from rotations around single bonds. 
However, in less exposed lysine side chains their rotational 
freedom can be restricted by the surroundings. In the solution 
structure of PpcA, for some of the lysine residues, like Lys^, many 
conformations were observed, while for Lys ' only a few 
conformations were observed (PDB codi;, 2LDO [1 1]). The 
distance between the weU-defined Cp of the selected lysine residues 
to the nearest iron(s) in the heme(s) is shown in Table 1 (see also 
Fig. 1). ^ 

To determine whether the difference in side chains can alter the 
local structure, the lysine residues were replaced in-silico by 
glutamic acid. The modeling studies show that the surface of 
lysines 9, 18, 22 and 60 can be replaced by glutamic acid (or 
glutamine) without any structural alterations. Even replacing the 
partially buried lysines 43 and 52 by glutamic acid did not indicate 
any steric clashes that cannot be avoided by slight adjustment of 
the standard side chain rotamers. Therefore, it can be concluded 
that the change observed in redox properties of mutants K43, K52 
and K60 (see below) are primarily due to replacement of the 
positive charge at that position. 



Screening the effect of the mutations on the relative 
heme oxidation fractions 

The heme oxidation fractions in the wild-type cytochrome were 
previously determined from the chemical shifts of the heme 
methyls 12'CH3', T'CHj™, 12'CH3'^ measured in die different 
oxidation stages from 2D-'H EXSY NMR experiments, according 
to the methodology summarized in the Materials and Methods 
section. Since PpcA exhibits fast intramolecular and slow 
intermolecular electron exchange on the NMR time scale, the 
individual heme signals in different oxidation stages can be 
discriminated and used to monitor the oxidation profiles of the 
hemes (for a review see [13,17]). The data obtained for PpcA 
showed that the heme oxidation fractions are modulated by the 
solution pH in the region 6—8. This modulation, known as redox- 
Bohr effect, is observed in the physiological pH range for G. 
sulfurreducens cellular growth giving the protein the necessary 
properties to couple e"/!!"*" transfer [17,20]. Thus, in the present 
work we first screened the oxidation profiles of the individual 
hemes in each mutant at pH 6 and 8 by measuring the chemical 
shifts of heme methyls 12'CH3\ 7'CH:/", 12'CH:/^ in the 
different oxidation stages. This is illustrated in the 2D-'H EXSY 
NMR spectra obtained at pH 8 for each mutant in Figure S5. 
Such strategy allowed us to identify and select the mutants with 
considerable differences in their heme oxidation profiles compared 
to the wild-type for further detailed thermodynamic characteriza- 
tion. 

The oxidation pattern of individual hemes of the mutants and 
wild-type proteins obtained at pH 6 and 8 are indicated in 
Figure 2 and 3. The comparison between the heme oxidation 
fractions, at both pH values, shows that the heme oxidation 
profiles were only slighdy altered for K9, K18 and K22 but 
significantiy changed for K43, K52 and K60 (c/. dashed and solid 
lines in Fig. 2 and 3). In the latter group of mutants (K43, K52 and 
K60) two main features emerge from the analysis of their oxidation 
patterns: (i) the hemes located near the replaced residue are the 
most affected ones (heme 111 for K60 and heme IV for K43 and 
K52 mutants) and (ii) the more affected hemes have higher 
oxidation fraction values compared to the wild-type protein {cf. 
dashed and solid hnes in Fig. 3). This latter observation can be 
rationalized on a pure electrostatic basis. Indeed, the replacement 
of a positively charged lysine with a neutral glutamine or with a 
negatively charged glutamic acid is expected to stabilize the 
oxidized form of the nearest heme by lowering its reduction 
potential and, concomitantiy, increase its oxidation fraction, 
compared to the other heme groups. Since K43, K52 and K60 
mutants showed significant differences in their oxidation profiles 
compared to the wild-type protein, their functional properties were 
characterized in detail in the present work. 

Redox characterization of heme groups and redox-Bohr 

centre in K43, K52 and K60 mutants 

In order to probe the functional role of each lysine residue in the 
modulation of the redox properties of PpcA, the detailed 
thermodynamic characterization of the six selected mutants 
(K43Q, K43E, K52Q, K52E, K60Q, and K60E) was carried 
out. The thermodynamic parameters of the wild-t)'pe protein were 
previously determined by fitting the pH dependence of the 
chemical shifts of heme methyls 12'CH3', 7'CH3^'' and 
12'CH3'^, measured in different stages of oxidation, together 
with data from visible redox titrations obtained at pH 7 and 8, 
within the framework of an electrostatic model that considers four 
interacting centres: three hemes and one protonatable centre [17]. 
In the present work, the same set of heme methyl groups were used 
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Table 1. Distance between the Cp 
code, 2LD0 [11]). 


of selected lysine residues to the nearest heme iron(s) in the solution structure of PpcA (PDB 
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to characterize the redox properties of PpcA mutants. 2D- H 
EXSY NMR spectra were collected at difierent pH values and the 
chemical shifts of heme methyls I2'CH3\ y'CH;/'' and 12'CH:/^ 
were measured for oxidation stages 1-3 (Fig. 4). As for the fully 
reduced form (Table SI), in the fully oxidized state (stage 3) the 
chemical shifts of the heme methyls from PpcA mutants are little 
affected (c/. corresponding black and red solid lines in Fig. 4). On 
the other hand, at intermediate oxidation stages (stages 1 and 2) a 
different scenario is observed. In this case, the chemical shifts of 
the heme methyls are quite difiFerent (Fig. 4). The most affected 
heme methyl in K43Q/E and K52Q,/E is 12'CH3'^, whereas in 
K60Q./E is y'CH:/". In each case, the most affected heme is 
located nearest to the mutated residue in the structure (see Fig. 1). 
By having essentially unaffected chemical shifts in the fuUy 
reduced and fully oxidized forms, the variations on the chemical 
shifts at intermediate oxidation stages reflect important changes in 
the heme oxidation fractions, i.e. changes in their redox 
parameters, during the oxidation of the proteins (see Fig. 4). It is 
important to note that in each oxidation stage, the number of 
oxidized hemes is fixed, i.e. one in oxidation stage 1 and two in 
oxidation stage 2 (Figure S2 - see Materials and Methods). Thus, if 
one particular heme is more oxidized, the others are necessarily 
less oxidized. This is clearly illustrated in the oxidation pattern of 
the heme groups at pH 6 and 8 indicated in Figure 2 and 3. 

In order to quantify the effect of the mutations on the redox 
properties of the heme groups, the thermodynamic model was fit 
to the pH dependence of the observed chemical shift of heme 
methyls, together with data from visible redox titrations (see 
Materials and Methods). The thermodynamic parameters deter- 
mined for the six mutants are indicated in Table 2. The quality of 
the fittings obtained for the pH dependence of the paramagnetic 
chemical shifts and for the visible redox titrations clearly shows 
that the experimental data is well described by the model (see 
black solid lines in Fig. 4). In the six mutants, the redox and redox- 
Bohr interactions are weaker compared to those of PpcA 
(Table 2). However, as observed for PpcA, the smallest redox 
interactions are established between the hemes that are structur- 
ally further apart (hemes I and IV) and the strongest redox-Bohr 
interactions with heme IV. The lower values of the redox 
interactions, in the context of a conserved heme core, might have 
their origin on structural rearrangements of charged groups that 
cause variations in local dielectric constants [28]. The positive 
values of the redox interactions indicates that the oxidation of a 
particular heme renders the oxidation of its neighbors more 
difiicult, which is expected on an electrostatic basis. In fact, the 
oxidation of one particular heme, i.e. transfer of one electron, 
decreases the number of negative charges in the vicinity of its 



neighbor, stabilizing the reduced form of the latter. Similarly, the 
negative redox-Bohr interactions (between the hemes and the 
redox-Bohr centre) indicate that the oxidation of the hemes 
facilitates the deprotonation of the redox-Bohr centre and vice- 
versa. The higher redox-Bohr interaction observed with heme IV 
in the mutants, indicates that the redox-Bohr centre remains 
associated with heme IV, as observed in the wild-type protein 
[10,1 1,17]. In relation to the heme reduction potentials, the most 
aflfected ones were those of hemes III (K60Q/E) and IV (K43Q,/E 
and K52Q/E). In all cases the reduction potential in the fully 
reduced and protonated proteins decreased, as expected from the 
removal of a positive charge in the mutant proteins from the 
neighborhood of the affected heme groups. The observed 
variations suggest an important role for residues Lys'''^, Lys''^, 
and Lys**" in the regulation of the redox properties of hemes III 
and IV. 

Impact of the mutations on the heme oxidation order at 
physiological pH 

As shown above, the reduction potential of each heme group is 
affected by the oxidation state of neighboring hemes (redox 
interactions) and by the solution pH (redox-Bohr interactions). To 
evaluate the effect of the mutations on the heme midpoint 
reduction potentials {eapp, i.e. the point at which the oxidized and 
reduced fractions of each heme are equally populated) at 
physiological pH, the oxidation curves of the individual hemes 
for K43, K52 and K60 mutants were computed from the 
thermodynamic parameters listed in Table 2 (Fig. 5 and Table 3). 
As for the wild-type protein, the redox interactions modulate the 
electron affinity of each heme so that the curves are non- 
Nernstian, showing some cross over during the redox titration 
(Fig. 5). The heme oxidation profiles of the mutants are different 
compared to the wild-type protein (c/. sohd and dashed lines in 
Fig. 5), which reflects the differences in the redox parameters and 
particularly in their heme reduction potentials. From the analysis 
of the individual oxidation profiles of PpcA mutants it is clear that 
at physiological pH, the hemes located in the proximity of the 
mutated residue are the most affected ones, i.e., heme III for K60 
and heme IV for K43 and K52 mutants. Compared to the wild- 
type the e^pp values are smaller, as expected from the replacement 
of a positive charge in their vicinity. A significant decrease in the 
reduction potential is observed for heme FV (34/54 mV and 49/ 
51 mV for K43Q,/E and K52Q,/E mutants, respectively). The 
significant decrease in the reduction potential of heme IV in K43E 
and K52Q/E mutants alters the heme oxidation order in these 
proteins (IV-I-III) compared to the wild-type (I-IV-III). The 
diflerence between the e„pp values of the first two hemes (heme IV 



PLOS ONE I www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | el 05566 



Role of Key Residues in Triheme Cytochrome PpcA from G. sulfurreducens 




Oxidation stage Oxidation stage Oxidation stage Oxidation stage 



Figure 2. Oxidation fraction of K9, K18 and K22 PpcA mutants (solid symbols and lines) and PpcA (open symbols and dashed lines) 

at pH 6 and pH 8. Data for Inemes I, III and IV are colored in green, orange and blue, respectively. The heme oxidation fractions were calculated 
according to equation x,- = (5i-5o)/(83-6o), where S,, Sq, and 63 are the observed chemical shifts of each methyl in stage i, 0, and 3, respectively. 
doi:1 0.1 371/journal.pone.01 05566.g002 



and I) is of the same order of magnitude compared to that 
observed for tlie wild-type protein (Table 3). However, tlie 
difference between the eapp values of the second and third heme 
to be oxidized is considerable higher (Table 3). On the other hand, 



the replacement of the positively charged lysine at position 60 
decreases significantly the eapp value of heme III (35 mV and 
38 mV for KBOQ and K60E, respectively) compared to the wild- 
type. Therefore, the oxidation curve of heme III is shifted to more 
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Oxidation stage Oxidation stage Oxidation stage Oxidation stage 



Figure 3. Oxidation fraction of K43, K52 and K60 PpcA mutants (solid symbols and lines) and PpcA (open symbols and dashed 
lines) at pH 6 and pH 8. Data for hemes I, III and IV are colored in green, orange and blue, respectively. The heme oxidation fractions were 
calculated as described in legend of Fig. 2. 
doi:1 0.1 371 /journal.pone.01 05566.g003 



negative values towards those of the other hemes (Fig. 5). As a 
consequence, in these proteins heme III is no longer the last one to 
be oxidized so that the heme oxidation order is I-III-IV for K60Q_ 
and (III,I)-IV for K60E. 



Role of K43, K52 and K60 on the functional mechanism of 
PpcA 

As shown above, the oxidation profile of the redox centres is 
highly dependent on the nature of the side chain at positions 43, 
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Figure 4. Fitting of the thermodynamic model to the experimental data for PpcA mutants. The blacl< solid lines are the result of the 
simultaneous fitting of the NMR and visible data. The three upper panels show the pH dependence of heme methyl chemical shifts at oxidation 
stages 1 (A), 2 (□), and 3 (O). The lower panel corresponds to the reduced fractions determined by visible spectroscopy at pH 7 (O) and pH 8 (□). 
The open symbols and the filled symbols represent the data points in the reductive and oxidative titrations, respectively. The solid red lines in each 
panel represent the best fit for the wild-type protein [17]. 
doi:1 0.1 371/journal.pone.01 05566.g004 



52 or 60. Thus, to evaluate the effect of each mutation on PpcA 
functional mechanism, the relative contribution of the 1 6 possible 
microstates (see Fig. S2) was determined as function of the solution 
potential for K43, K52 and K60 mutants (Fig. 6). Such studies 
were previously undertaken for the wild-type protein and a 
coherent electron transfer pathway coupled to proton transfer was 
established [17]. The relative variation of the microstates in PpcA 
is also depicted in Figure 6 together with those obtained in the 
present work for K43, K52 and K60 mutants. While in the wUd- 
type protein a concerted e /H"^ transfer occurs between oxidation 
stages 1 and 2, it is clear that the relative contributions of different 
microstates to the overall population in each stage are significantly 
altered in the lysine mutants (Fig. 6). In fact, in the wild-type 
cytochrome the oxidation stage 0 is dominated by the fuUy 
reduced and protonated form Pofj and stage 1 is dominated by the 
oxidation of heme I (Pih) while keeping the redox-Bohr centre 
protonated. Oxidation stage 2 is dominated by the oxidation of 
hemes I and IV and deprotonation of the acid-base centre {P14), 
that remains deprotonated in stage 3 upon full oxidation of heme 
III {P134). Therefore, a route is defined for electron transfer in 
PpcA: Pq^ Pih ^ P14 Pi 34- This is not the case for the 
K60 mutants for whom the lowering of the eapp of heme III brings 
the midpoint reduction potential values of all hemes groups closer 
to each other and favours the oxidation of this heme at lower 
reduction potentials (c/. dashed and sohd hnes on Fig. 5). 
Consequently, for the K60 mutants, two microstates {Pih and 
P3//) are now dominant in the first oxidation stage and, thus, no 
preferential pathway for electron transfer is observed (Fig. 6). On 
the other hand, for K43 and K52 mutants a different scenario is 
observed. As mentioned above, the major effect of replacing the 
positive charge is the lowering of the eapp values of heme IV. 
Consequentiy, this heme is more oxidized in the first oxidation 
stage, compared to PpcA. In the case of K43E, K52Q,and K52E 
the significant reduction on the e^pp of heme IV alters the 



oxidation order of the hemes (see Fig. 6 and Table 3) so that heme 
IV is now the first heme to oxidize. In all these mutants heme III is 
clearly the last one to oxidize. Therefore, and since the eapp of 
heme I is closer to that of heme IV (cf. solid green and blue lines in 
Fig. 5) the microstate with heme I and IV oxidized {P14) in the 
mutants has a much higher population compared to the wild-type 
protein (Fig. 6). In the mutants, the oxidation stage 0 is also 
dominated by the protonated form Pqh, but the microstates of 
oxidation stage 1 are overcome by the P14 curve, which intercepts 
earlier the curve of microstate Poh- The microstate P24 dominates 
the oxidation stage 2, whereas Pi 34 dominates stage 3. Thus, in 
these three mutants, a different preferential route for electrons is 
established, favoring a proton-coupled 2e~ transfer step between 
oxidation stages 0 and 2: Pqh P14 P134 as observed for 
PpcD and M58K and M58N mutants [17,29]. 

Taken these observations together, it is clear that residues Lys*'^, 
Lys' and Lys ' modulate the reduction potential of their closest 
heme, which in turn controls the microscopic redox states that can 
be accessed during the redox cycle of the proteins. However, the 
repercussion of the substitution of a positive charge close to heme 
III (Lys''° mutants) or IV (Lys*'' and Lys''^ mutants) has distinct 
effects on the functional mechanism of the protein. In case of 
Lys''", the replacement of the positive side chain by the neutral or 
negatively charged one altered the reduction potential of heme III 
in a way that it is no longer the last heme to oxidize. In this case, 
no preferential pathway for e /Yt' transfer is observed. In 
contrast, in the case of Lys*'^ and Lys^^ mutants, the removal of a 
positively charged side chain either at position 43 or 52 lead to a 
higher stabilization of the oxidized form of heme IV, compared to 
the other two hemes, and a different functional mechanism 
emerges in which a 2e step coupled with proton transfer is 
observed. This suggests that the reduction potential of heme III 
relative to the other two hemes seems to be crucial in enabling 
these proteins to couple electron transfer with deprotonation of the 
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Table 2. Thermodynamic parameters for PpcA mutants. 



Energy (meV) 



Protein 


Heme 1 


Heme III 


Heme IV 


Redox-Bohr centre 


K43Q 


Heme 1 


-161 (5) 


20 (2) 


6(3) 


-17(4) 


Heme III 




146 (4) 


30 (3) 


-15 (4) 


Heme IV 






141(5) 


-42 (4) 


Redox-Bohr centre 








463 (8) 


K43E 


Heme 1 


161 (3) 


16 (2) 


-3 (2) 


-11 (3) 


Heme III 




145 (3) 


20 (2) 


-6 (3) 


Heme IV 






167 (3) 


-38 (3) 


Redox-Bohr centre 








463 (6) 


K52Q 


Heme 1 


154 (3) 


20 (2) 


0(2) 


-15 (3) 


Heme III 




138 (3) 


20 (2) 


-11 (3) 


Heme IV 






160 (4) 


-43 (3) 


Redox-Bohr centre 








466 (6) 


K52E 


Heme 1 


154 (3) 


20 (2) 


-1 (2) 


-14(3) 


Heme III 




142 (3) 


27 (2) 


-22 (3) 


Heme IV 






160 (4) 


-59 (3) 


Redox-Bohr centre 








482 (6) 


K60Q 


Heme 1 


164 (3) 


22 (2) 


6(2) 


-15 (3) 


Heme III 




162 (3) 


30 (2) 


-16 (3) 


Heme IV 






-132(3) 


-40 (3) 


Redox-Bohr centre 








461 (6) 


K60E 


Heme 1 


-151 (5) 


24 (3) 


9(3) 


-23 (5) 


Heme III 




-157 (6) 


35 (3) 


-29 (5) 


Heme IV 






121 (6) 


-49 (5) 


Redox-Bohr centre 








482 (11) 


PpcA 


Heme 1 


154 (5) 


27 (2) 


16(3) 


-32 (4) 


Heme III 




138 (5) 


41 (3) 


-31 (4) 


Heme IV 






-125 (5) 


-58 (4) 


Redox-Bohr centre 








495 (8) 



For comparison the values previously obtained for PpcA [17] were also included. For each cytochrome, the fully reduced and protonated protein was taken as reference. 
Diagonal values (in bold) correspond to oxidation energies of the hemes and deprotonating energy of the redox-Bohr centre. Off-diagonal values are the redox (heme- 
heme) and redox-Bohr (heme-proton) interactions energies. All energies are reported in meV, with standard errors given in parenthesis. 
doi:l 0.1 371 /journal.pone.OI 05566.t002 



redox-Bohr centre. Indeed, only when heme III is the last one to 
oxidize (higher eapp value) preferential e /H"*" pathways are 
established, which are reinforced by a higher separation between 
the eapp values of the second and third hemes to oxidize (see 
Table 3). However, the pathway varies with the separation 
between the e^pp values of heme III and its predecessor in the 
order of oxidation (Table 3). In the case of PpcA such separation 
was 18 mV and the route for electron transfer was: Pqh Pjh^ 
P]4 Pi34- In case of K43Q_the same route was observed but the 
slightly higher separation between the eapp values (33 mV versus 
1 8 mV in the wild-type) increased the contribution of microstate 



P]4. Finally, in K43E, K52Q, and K52E the higher separation 
between the eapp values of hemes III and I (48, 46 and 40 mV, 
respectively) lead to a significant contribution of the microstate Pm 
so that a drEFerent preferred route for electrons was observed: Pqh 
Pi4 Pi34- Interestingly, the results obtained from the 
thermodynamic characterization of the G. sulfurreducens PpcA 
family members (PpcA, PpcB, PpcD, and PpcE) [17], as well as 
PpcA mutants of residue Met' , a residue that controls the solvent 
accessibility of heme III [29], have shown that preferential 
pathways for e"/!!""" coupling were established only when heme 
III has the highest reduction potential (Table 3). 
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Figure 5. Oxidized fractions of the individual Kiemes for PpcA mutants (solid lines) and wild-type (dashed lines) at pH 7.5. The curves 
were calculated as a function of the solution reduction potential using the parameters listed in Table 2. The midpoint reduction potentials of the 
hemes (e^pp) are also indicated. 
doi:1 0.1 371/journal.pone.01 05566.g005 
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Figure 6. Molar fractions of the 16 individual microstates (described in Fig. S2) of PpcA mutants and the wild-type cytochrome at 
pH 7.5. The curves were calculated as a function of the solution reduction potential using the parameters listed in Table 2. Solid and dashed lines 
indicate the protonated and deprotonated microstates, respectively. For clarity only the relevant microstates for a given mutant are labelled. 
doi:10.1371/journal.pone.0105566.g006 



Finally, the impact of the changes in the heme redox potentials, 
heme-heme and redox-Bohr interactions on the macroscopic 
behaviour of K43, K52 and K60 mutants were evaluated. The 
global oxidation curves of these mutants are depicted in Fig. 7. In 
all cytochromes, the E^pp values [i.e. the point at which the 
oxidized and reduced fractions of each protein are equal) are more 
negative compared to PpcA indicating that the functional working 
potential ranges in the mutants are shifted to lower redox potential 
ranges. The lower E„pp values of the mutants are positioned to 
thermodynamically favour the reduction of downstream redox 
partners in the Gs respiratory chain. 



Conclusions 

In the present work we studied the roles of six positively charged 
lysine residues on the redox properties of the triheme cytochrome 
PpcA. The selected residues are located in the vicinity of the heme 
groups and were replaced separately with glutamine and glutamic 
acid to evaluate the eflFects of a neutral and negatively charged 
residue at each position. The results showed that Lys (near heme 
IV), Lys'" (near heme I) and Lys^^ (between hemes I and III) have 
no effect on the redox properties of the heme groups and are 
probably involved redox partner recognition. However, this 
hypothesis can only be fuUy confirmed upon identification of 
PpcA redox partners. On the other hand, for Lys and Lys' 
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Figure 7. Macroscopic oxidation curves of PpcA mutants (black solid lines) and wild-type (red dashed line) at pH 7.5. The curves were 
calculated as a function of the solution reduction potential using the parameters listed in Table 2. The midpoint reduction potentials of proteins 
(fapp) are indicated in the inset. 
doi:1 0.1 371/journal.pone.01 05566.g007 
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(both near heme TV) and Lys (near heme III) the results showed 
that these residues are crucial in the selection of microstates that 
allow the protein to establish preferential e /H^ transfer pathways 
and, therefore, in the regulation of the functional mechanism of 
PpcA. Indeed, the removal of positive charges at positions 43, 52 
or 60 altered the redox properties of the hemes and the balance of 
the global network of cooperativities. The data obtained suggests 
that the preferred e /IT^ transfer pathway observed for PpcA is 
strongly dependent on the reduction potential of heme III relative 
to the other hemes. In fact, when the reduction potential of heme 
III is similar to the reduction potentials of the other two hemes, as 
in the case of K60 mutants, the preferential pathway for electron 
transfer is disrupted. In contrast, when heme III is clearly the last 
one to oxidize in native, K43 or K52 mutants, there is a well- 
established e /H"*" transfer pathway that involves a single or a two 
electron transfer coupled to proton transfer. The latter is favored 
by a larger diflFerence between the reduction potential of heme III 
and its predecessor in the order of oxidation. Taken together the 
information obtained in the present and previous studies, we show 
that this feature is independent from the order of oxidation of the 
first two hemes and that the functional mechanism of PpcA relies 
on a fine tuned balance of redox and redox-Bohr interactions to 
assure a coherent electron transfer pathway txjupled to proton 
transfer that allows the protein to perform e~/H^ energy 
transduction. The detailed characterization of the redox properties 
of PpcA mutants wiU enable rational design of functional PpcA- 
based proteins operating at different working potential windows, 
but, more importantly, preserving the coherent electron transfer 
pathways that contribute to cellular energ)^ transduction. Having 
been successful in the design of proteins that work at more 
negative working potential ranges, with the concomitant increase 
of the driving force for electron transfer, this work sets the 
foundations to develop strains carrying mutant cytochromes 
rationally designed to increase respiratory capabilities of Geohacter 
species. 

Supporting Information 

Figure SI Alignment of G. sulfurreducens triheme 
cytochromes PpcA-E amino acid sequences. The conserved 
residues are boxed: non-heme attached in black and heme 
attached in gray, with the corresponding heme number shown 
below in Roman numerals. For cytochromes PpcB-E, the 
percentage of sequence identity in relation to PpcA is indicated. 
The percentage of lysine residues in each protein are shown in 
parentheses. 
(TIF) 

Figure S2 Electronic distribution scheme for a triheme 
cytochrome with a proton-linked equilibrium showing 
the 16 possible microstates. The green and orange circles 
correspond to the protonated and deprotonated microstates, 

respectively. Inner circles represent heme groups, which can be 
either reduced or oxidized and are colored gray or white, 
respectively. The microstates are grouped, according to the 
number of oxidized hemes, in four oxidation stages connected 
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